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ABSTRACT: The objective of this study was to get a comprehensive understanding of how genes in chicken shell gland modulate 
eggshell strength at the early stage of active calcification. Four 32-week old of purebred Xianju hens with consistent high or low shell 
breakage strength were grouped into two pairs. Using Affymetrix Chicken Array, a whole-transcriptome analysis was performed on 
hen's shell gland at 9 h post oviposition. Gene ontology enrichment analysis for differentially expressed (DE) transcripts was performed 
using the web-based GOEAST, and the validation of DE-transcripts was tested by qRT-PCR. 1,195 DE-transcripts, corresponding to 
941 unique genes were identified in hens with strong eggshell compared to weak shell hens. According to gene ontology annotations, 
there are 77 DE-transcripts encoding ion transporters and secreted extracellular matrix proteins, and at least 26 DE-transcripts related to 
carbohydrate metabolism or post-translation glycosylation modification; furthermore, there are 88 signaling DE-transcripts. GO term 
enrichment analysis suggests that some DE-transcripts mediate reproductive hormones or neurotransmitters to affect eggshell quality 
through a complex suite of biophysical processes. These results reveal some candidate genes involved with eggshell strength at the early 
stage of active calcification which may facilitate our imderstanding of regulating mechanisms of eggshell quaUty. (Key Words: 
Chicken, Microarray, Differentially Expressed Genes, Eggshell Strength) 



INTRODUCTION 

The chicken eggshell is a porous bioceramic container 
which protects the egg against physical damage and 
microbial contamination. Avian eggshell consists of the 
innermost bilayered membranes, a calcified layer composed 
of a mamillary and pallisade layer, and the outermost 
cuticle. The calcified layer consists of both inorganic 
minerals and extracellular matrix. It is well known that the 
shell mineral amount (thickness or density) is the main 
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factor contributing to the mechanical properties of the 
eggshell (Ahmed et al., 2005). However, the organic matrix, 
although its content in the calcified layer is only 2 to 3.5%, 
is of great importance to the deposition of bicarbonate and 
calcium ions, and to eggshell strength by controlling calcite 
crystal nucleation, growth, size and orientation (Greenfield 
et al., 1984). 

The organic matrix in the calcified layer is comprised of 
a complex suite of components. In the acid soluble part of 
chicken eggshell matrix, 520 proteins have been identified 
(Mann et al., 2006), including several abundant proteins 
such as ovalbumin (Hincke, 1995), ovotransferrin (Gautron 
et al., 2001b), lysozyme (Hincke et al., 2000), osteopontin 
(Pines et al., 1995), sialoprotein (Solomon, 1999), clusterin 
(Mann et al., 2003), ovocleidin-17 (Hincke et al., 1995), 
ovocleidin-23 (Mann, 1999), ovocleidin-1 16 (Carrino et al., 
1997), ovocalyxin-32 (Gautron et al., 2001a) and 
ovocalyxin-36 (Gautron et al., 2007). Many of the above 
components have been reported to undergo various post- 
translation modifications, which allow them to be effective 
chelators for interacting with the inorganic materials (Veis, 
1989; Reyes-Grajeda et al., 2004; Mann et al., 2007), or to 
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mediate protein-protein interactions to facilitate the 
assembly of the organic matrix (Lakshminarayanan et al., 
2002; Ney et al., 2006). 

It has been demonstrated that some genes in hen oviduct 
are associated with eggshell formation, whose expression is 
dependent on mechanical strain (Pines et al., 1995; Lavelin 
et al., 1998; Lavelin et al., 2002). It is proposed that some 
genes may function as crucial modulators for eggshell 
quality through regulating signal transduction, ion 
transportation, expression or modification of organic 
components, and many other processes. However, despite 
the importance of eggshell strength in the poultry industry, 
very few transcriptome-wide studies regarding this trait 
have been published to date (Yang et al., 2007; Dunn et al., 
2009; Jonchere et al., 2010). 

It is well documented that various parts of the avian 
eggshell are formed in specific regions of the oviduct as the 
egg passes through them. During the laying sequence, about 
4 h after previous oviposition, the next egg arrives at and 
will take about Ih to pass through the white isthmus, in 
which the bilayered shell membranes are built around the 
egg. Then the egg enters the initial part of the shell gland, 
the red isthmus (tubular shell gland), and stays there for 
about 5 h to form mammillary knobs (Reyes-Grajeda et al., 
2004). Finally the egg reaches the uterus (the main part of 
shell gland) and stays for an additional 15 h to form the 
palisade layer (Creger et al., 1976). It is known that the 
mamillary layer is the base of calcite crystal nucleation and 
crystal growth, and the palisade layer is the main part of the 
calcified shell, both of which affect global eggshell quality 
(Reyes-Grajeda et al., 2004; Jonchere et al., 2010). 

In this study, we focused on the shell gland (uterus 
tissue near the red isthmus) at about 9 h post oviposition 
(corresponding to the early stage of active calcification, or 
to the transition stage from mammillary knob formation to 
construction of the palisade layer), and identified 
differentially expressed genes (DE-genes) in the layers with 
high shell strength compared to those with weak eggshell. 
Our results provide insight into the candidate genes 
involved in the mamillary layer formation and calcification 
that is crucial to the mechanical properties of avian 
eggshells. 

MATERIALS AND METHODS 

Animal treatments 

Ninety purebred Xianju hens (a widely-bred Chinese 
indigenous chicken breed) of 28 weeks old were 
individually housed in laying cages. Birds were maintained 
under a cycle of 16 h light and 8 h dark. All birds were fed 
ad libitum with water and a mash layer diet (165 g protein, 
35 g Ca, 11.29 MJ ME/kg, as recommended by NRC of 
China, 2004). 



After 10 d of adaptation for hens, the oviposition time 
of each egg was initiated to be observed and recorded, then 
egg weight and shape index (length/width) were measured 
immediately. Following strength testing, the egg content 
was discarded and the shell was washed, dried at room 
temperature and weighed. Shell thickness without 
membranes was measured with a digital micrometer. Shell 
index (g/100 cm^) (Sauveur, 1988) was calculated as I = 
(C/S)xlOO, in which C is the weight of shell with 
membranes, S is the shell surface (cm^) with S = 4.68xP^'^ 
where P = egg weight (g). All above measurements were 
consecutively carried out daily for 16 d. 

Finally, 2 groups of 2 hens with consistent high or low 
shell breakage strength were found. The differences 
between the eggshell properties of the selected 4 hens were 
analyzed by One-way ANOVA variance analysis in SPSS 
statistic software. 

The four hens of interest were humanely sacrificed 
about 9 h after the previous oviposition. It is of note that all 
of the sacrificed hens had eggs in their uteruses (Figure 1 A). 
The fat was removed from the uterus tissues near red 
isthmus and the tissues were then frozen in liquid nitrogen 
immediately and stored at -80°C. The animal treatments 
were approved by the Commission for Animal Welfare of 
Zhejiang A&F University. 



in IV 



Figure 1. Eggs and forming eggshells obtained from hen shell 
glands at the moment of tissue sampling. Figure lA: Eggs 
obtained from the shell glands of hens when sampling utems 
tissues. A-I is the egg from #35 hen, A-II from #19 hen, A-III from 
#80 hen, A-IV from #40 hen; these hens were all slaughtered at 
about 9 h after previous oviposition (P.O.). While A-V and A- VI 
eggs are from another 2 hens culled at 11.5 to 12 h after P.O., 
respectively. Figure IB: The forming "shell" sampled at different 
stage of eggshell formation. B-I is the "shell" from the hen 
slaughtered at 4 h after P.O.; B-II is the "shell" of above egg from 
#40 hen. B-III is the "shell" sampled at 12 h post oviposition. 
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Measurement of eggshell strength 

After egg weight and shape index measurements, the 
uncracked fresh eggs were individually placed lengthways 
with its blunt end upward in the FHK testing machine 
(Fujihara Co., Tokyo, Japan), and the vertical pressure was 
increasingly loaded upon the eggshell until the eggshell 
cracked and the eggshell strength was recorded as the 
maximum load (kgf). 

RNA preparation 

About 500 mg of the tissue of the uterus near red 
isthmus, including the mucosa, muscularis and outer serosa 
was powdered under liquid nitrogen. The total RNA was 
extracted using the RNAiso Plus Mini Kit (TaKaRa, Dahan, 
P.R. China) according to the manufacturer's instructions. 
RNA concentration and purity were measured by a 
NanoDrop spectrophotometer (NanoDrop® ND-1000, 
NanoDrop Technologies, DE). 

Microarray hybridization and image acquisition 

Microarray analysis was performed by the Bioassay 
Laboratory of CapitalBio Corporation (CapitalBio Co., 
Beijing, China). Briefly, the RNA integrity was firstly 
assessed using a Bioanalyzer (Agilent Technologies, 
Cheshire, UK), then 2 )ig of total RNA was used for reverse 
transcription and biotin-labeled cRNA synthesis according 
to the manufactures' instructions, and finally subjected to 
microarray hybridization. The Affymetrix GeneChip® 
Chicken Genome Array (Affymetrix, Santa Clara, CA, 
USA) was used in this study, which contains 38,535 
probesets corresponding to >28,000 chicken genes. 
Following 16 h of hybridization, the arrays were 
immediately washed, stained and scanned using 
Affymetrix® GeneChip® scanner 3000 (Affymetrix, Santa 
Clara, CA, USA), and the image files were processed into 
raw CEL intensity files using GeneChip Operating Software 
(GCOS version 1.2). 

Pre-processing and normalization of microarray data 

The raw intensity files generated by GCOS were 
imported and processed by R with Bioconductor packages. 
The total RNA quality was firstly verified statistically again 
by plotting the 5 '-3' hybridization signal trends across all 
target transcripts. Then the microarray intensity was 
processed into transcript expression by the Affymetrix 
MAS5.0 method implemented in the R package, a 
procedure including background normalization, PM/MM 
probe correction, expression summarization and constant 
normaUzation on probeset level. 

Identification of DE-transcripts 

To identify DE-transcripts, the 4 array samples were 
first grouped into two pairs of high vs. low eggshell 



strength according to eggshell property differences of the 
hens (see results). According to Cheuk and Cheng (2011), 
Affymetrix platform is relatively precise and sensitive in 
detecting signals, the DE-transcripts were identified as 
those with fold-change >= 2 in either of the two pairs of 
comparison and a statistical significant difference between 
high strength and low strength samples (p<0.05, Welch t- 
test). It is of note that the log-odds values (Lods) of 
expression fold-change were used in the analysis; therefore, 
the DE-transcripts always have an absolute Lods value no 

less than 1 ( |LoJ5| > 1 ). 

Gene ontology enrichment analysis 

Gene ontology enrichment analysis for DE-transcripts 
was performed using the web-based GOEAST (Zheng and 
Wang, 2008) Affymetrix analysis tool, with FDR cut-off of 
0.05 using YekutieU's FDR adjustment method. 

VaUdation of differential expression by qRT-PCR 
experiments 

Twenty-one DE-transcripts, with a fold-change ranging 
from low to high, were selected for further validation with 
qRT-PCR experiments; and all two groups of microarray 
samples were tested. Glyceraldehyde 3-phosphate 
dehydrogenase (GAPDH) was used as the internal reference 
in all the PCR experiments. The primer sequences for qRT- 
PCR experiments can be found in Table 1. 

To begin, total RNA was individually reverse 
transcribed with the SYBR® PrimeScript™ RT-PCR kit II 
(TaKaRa, Dahan, China) according to the manufacturer's 
instructions. Then above RT-PCR kit was further used for 
fluorescence detection on an ABI Prism 7500 Sequence 
Detection System (Applied Biosystems, USA). All samples 
were analyzed in triplicates. 

Dissociation curve analysis was conducted to ensure 
that a single PCR product with appropriate size was 
amplified in each reaction. On the other hand, the 
examination of PCR efficiency was performed based on 
LinRegPCR program (12.X) (Ramakers et al., 2003; Ruijter 
et al., 2009) to ensure internal and target transcript primers 
were amplified with similar efficiency. 

The differential expression levels (Log2 units) were 
calculated using the equation Log2 units (high versus low) = 
-AACt, where AACt = (Ct*-Ct''')-(Ct"-Ct''). Ct is the 
threshold cycle number when the amount of amplified 
product reaches a stable threshold. Ct* and Ct''' represented 
the Ct of target transcript and internal reference transcript of 
"high eggshell strength sample", respectively. 
Correspondingly, Ct" and Ct'' represented the Ct of target 
transcript and inner-reference transcript of "low eggshell 
strength sample", respectively. 
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Table 1. Descriptions of specific primers used for real-time RT-PCR 



Gene symbol 


Accession no. 


Forward primer (5'-3') 


Reverse primer (5 '-3') 


Amplicon (bp) 


ACYP2 


XM_ 


.419292 


CGGCTCGCTCAAGTCGGTGG 


GGCCCTGAACTTGGCCCGTC 


152 


AMDHDl 


XM. 


.416158 


GCACTGGGAAGTGCGTATTGCCA 


TCTTCCGTGGCCTTCCTGGTGT 


175 


ATP6V1A 


NM_ 


.204974 


TGCAACATGGCAGGTGCTGCT 


TGCCAGGCCCCAGTTCCACT 


187 


CA5B 


XM_ 


.414195 


CAGCTTGGCCACCTGCACTCG 


ACACGTCGCTGGGTCGTAGCT 


175 


CHST3 


NM. 


.205121 


TGATGGCCACCACACGCACC 


CTGCAGCACGTCGCGGTACA 


170 


CULIZAI 


NM. 


.205021 


AOOCGAG 1 CI 1 CCCCGACGG 


UCUC 1 G 1 CCl C- Al G 1 CI GCCC 


171 


A Tin 1 

CRABPl 


NM. 


.001030539 


CGCCCCGCCATGCCTAACTT 


AACTGGTCCCCGTCCTGGCG 


161 


CRYBBl 


NM_ 


.204180 


ACCTGGCGGACTGCGGGTT 


CGGTAGCTGCTGGACCAGGTG 


151 


tXULoB 


XM_ 


.420892 


AACCCCACCACAGCCCTCGT 


TGGCTGTTGATGAGGCCGCG 


149 


FGB 


XM. 


.420369.2 


GCTGCTCCTGCTGCTCCTGC 


GTGCCACGGGCCTGAGTGTG 


155 


GAS2L3 


XM_ 


.416172 


GGAGTAGTGCTGGCAGTCCTGC 


CCTGGGCCGTGTCTGGGAGT 


193 


GIT2 


NM. 


.204206 


TCGCTTGCCATGCCGTGAGG 


GCAACGTGGAGCGGGGTGTT 


168 


MAN1A2 


XM_ 


.416490 


ACGTGGACACCAGCAAGGGGG 


TCCTTTGCCTCTTCCAGGGCCTTT 


148 


NDST4 


XM_ 


.420638 


CGAGCAGCTTCCCTCATCCCCAA 


TGCCCAGGGGCTTGACGTAA 


156 


NPY 


NM. 


.204587 


GAGGACGCTCCCGCAGAGGA 


TCGAAGGGTCTTCAAACCGGGA 


175 


OC416916 


XM_ 


.415207 


TGGAGGTGGAGCACAAACATCTGC CCACCGAGCACACAGCCAGAAA 


200 


PLCXDl 


NM_ 


.001128637 


CCTGGCCTGCAGGAATTTTGATGG 


AGCCACGCTGCCACATGGTC 


137 


RCJMB04_34k20 


NM_ 


.001031112 


GGACAGGCGGGCGAGAGAGT 


TGGTGGTAACACGCACGCTGA 


126 


SLC8A1 


XM_ 


.415002 


CGTGTTTGTGGCACTGGGGACA 


ATGGCCGCGATGGACCAAGC 


159 


TBXASl 


XM. 


.416334 


TGTGTGGTGCTGGGACAGCGT 


ATACAGCCACGGGGTCCTGCT 


188 


WDR72 


XM_ 


.425069 


GGCTGTTATCAGGGGGCCAGGA 


GCACACGCAGCACACTACGC 


161 


GAPDH 


NM_ 


.204305 


GGGCTGCTAAGGCTGTGGGG 


TCAGGGGCCCATCAGCAGCA 


177 



RESULTS 

Eggshell quality of hens under study 

Among 90 tested hens, only 6 of them laid eggs at a 
similar laying rate with consistent high eggshell strength 
(defined as >4.5 kgf) or low strength (defined as <3.5 kgf). 
Two of these 6 hens were sacrificed at about 11.5 to 12 h 
post oviposition, but the eggshells collected from the shell 
glands showed more calcification extent than expected 
(Figure lA and IB). To focus on the initial stage of active 
calcification, we decided to use the uterus tissues near red 
isthmus from the remaining 4 hens, namely #19, #35, #40 
and #80, at about 9 hr post oviposition. The eggs and 
eggshells harvested from these 4 hens at the moment of 
tissue sampling are shown in Figure lA and Figure IB. The 
eggshell breaking strength is also shown in Table 2 for 
these 4 hens, with consistent high (#19 and #40) or low 
(#35 and #80) eggshell quality. 

To eggs from #19 and #35 hens, the differences of both 



shell strength and shell weight were very significant 
(p<0.01, t-test), but there was no significant difference 
(p>0.05) for other eggshell quality metrics, such as shell 
thickness or shell index (Table 2). On the other hand, the 
differences of all of above eggshell metrics between eggs 
from #40 and #80 hens were very significant (p<0.01) 
(Table 2). To get rigorous microarray data, therefore, we 
grouped the #19 with #35 hens due to the similarity of some 
of eggshell mechanical properties of the paired individuals; 
while #40 and #80 hens were also grouped as another pair. 

Differentially expressed transcripts 

The expression level of all probesets in 4 array samples 
were analyzed, and 1,195 DE-transcripts between uterus 
samples with high shell strength and low shell strength were 
identified. These DE-transcripts correspond to 941 unique 
genes. Among them, 407 genes were up-regulated in high 
strength samples comparing to low strength samples, and 
the other 534 genes were down-regulated. The expression 



Table 2. Parameters related to eggshell quality of hens in this study 


Hen shell strength (kgf) 


shell thickness (mm) 


shell index (g/100 cm^) 


shell weight (g) 


egg weight (g) 


shape index 


#40 5.17±0.40'^ 


0.367±0.016* 


8.07±0.28 ^ 


4.297+0.144" 


38.4611.52 


1.27410.024'' 


#19 4.75±0.21 


0.32810.012"^ 


7.1810.25" 


3.96810.195'' 


40.7211.26 ^ 


1.32210.030^ 


#35 2.99±0.71 ^ 


0.324±0.023 


6.9710.68 ^ 


4.36310.426 " 


48.6011.38 ^ 


1.27810.032'^'' 


#80 2.54+0.69 


0.272+0.022 


6.25+0.7 1 


3.45110.417 ' 


40.5811.06 


1.31310.052-^'' 



Values are from eggs laid by each hen of interest during the period of observation. Distinct capital letters in the same column indicate parameters between 
hens with a significant difference (p< 0.01), and distinct small letters indicate the significant difference is at level p<0.05. 



Liu et al. (2013) Asian Australas. J.Anim. Sci. 26:609-624 



613 



D 

0 5 10 

Value 




■>!j- T- 00 00 



Figure 2. Heatmap and dendrogram of differentially expressed 
transcripts (DE-transcripts). Each cell represents the (normalized) 
gene expression value for given DE-transcript (row) in the 
specified sample (column). Cell colors indicate gene expression 
level: red: highly expressed; yellow: medium expression; green: 
lowly expressed. Row-side and column-side dendrogram represent 
the hierarchical clustering of DE-transcript expression for 
different transcripts or samples, respectively. Clustering is based 
on "complete-linkage" method using Euclidean-distance. 

profile of all 1,195 DE-transcripts is shown in the heatmap 
in Figure 2. As shown in the heatmap, samples #19 with 
#40 and #35 with #80 were grouped as clusters among 
different samples, consistent with the similarity of their 
eggshell quality. 



According to gene ontology annotations, the DE- 
transcripts are involved in a variety of biological processes. 
The most prominent DE-transcripts were found related to 
the following processes: signal transduction (88 DE- 
transcripts), ion transport and extracellular matrix 
organization (77 DE-transcripts), carbohydrate metabolism 
and protein modification (26 DE-transcripts) (Table 3). 

Furthermore, avian calcified eggshell is a biomaterial 
composed of calcium salt and special ECM. The ECM is 
mainly comprised of coUagens, glycoproteins and 
proteoglycans. Among the DE-transcripts, COL8A2, 
C0L12A1, C0L13A1, LOC424798, LAMA2, LAMA4, 
LAMB4, and LAMCl may be related to extracellular matrix 
formation; while CHST3, GALNTLl, NDST4, LARGE, 
P0FUT2, RCJMB04_28123, and MAN1A2 are all localized 
in the endoplasmic reticulum or Golgi apparatus, and likely 
mediate the processes of carbohydrate metabolism, or 
posttranslation glycosylation modification. 

Gene ontology (GO) term enrichment of DE-transcripts 

It is of note that although many DE-transcripts were 
found related to various biological processes according to 
their ontology annotations, they are not necessarily 
coiTelated to the eggshell quality, due to random noise or 
other non-specific confounding factors commonly existing 
in microaiTay or other high-throughput experiments. 
Therefore, using web-based GOEAST (Zheng and Wang, 
2008) we further identified significantly enriched GO terms 
among all the DE-transcripts. According to biology 
processes or molecular functions, the enriched GO terms 
can be roughly classified into several groups (Tables 4 and 
5). 

A group of processes are involved in reproductive 
hormone regulation, which contain Somatotropin secreting 
cell differentiation (GO:0060126), adenohypophysis 
development (GO:0021984), and response to estradiol 
stimulus (GO:0032355) (Table 4). 



Table 3. i) DE-transcripts related with signaling, ion transportation, extracellular matrix protein, and carbohydrate metabolism or post- 
translation glycosylation modification 



Gene symbol 


Transcript ID 


Log2 units (strong VS weak) 


p-value 


Category 


LOC771699 


XM_001234946 


3.853 


0.0004 


signaling 


SH3PXD2A 


XM_421741 


3.410 


0.0397 


signaling 


LOC429955 


XM_427511 


3.296 


0.0034 


signaling 


PDCL2 


XM_420702 


3.100 


0.0115 


signaling 


LOC430487 


XM_428042 


3.055 


0.0074 


signaling 


SEMA3G 


XM_414289 


2.933 


0.0213 


signaling 


RH0BTB2 


XM_001232709 


2.905 


0.0069 


signaling 


RXFPl 


XM_420385 


2.795 


0.0205 


signaling 


PIK3C2B 


XM_417956 


2.777 


0.0194 


signaling 


OR10A7 


XM_425093 


2.768 


0.0067 


signaling 


NPY 


NM_205473 


2.401 


0.0329 


signaling 


PDE8B 


XM_425218 


2.401 


0.0498 


signaling 


GREM2 


XM_419552 


2.284 


0.0423 


signaling 
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Table 3. ii) DE-transcripts related with signaling, ion 
translation glycosylation modification (Continued) 


transportation, extracellular matrix protein, and carbohydrate metabolism or post- 


VjcllC ayiiiuui 


lldilsLIipL ILJ 


LjUgZ Uliila V^allUIlg V o WCoKJ 


1 JJ-VaiUC 




QUA /~^Ci 


AJVl_4zUUyo 


0 0 1 < 


n n/io'z 
U.U4zj 


signaling 


LUL-jyOjOJ 


JNJVI_zUj4UU 


O 1 


A nAAO 


signaling 


MArKBrl 


AK_Uzo / / z 


O 1 C7 
Z. 1 J / 


A ni on 
U.UliS/ 


signaling 


UaIK 


IN M_UU 1 1 joy 


z.izo 


A ni 'XA 

U.Ui j4 


signaling 


Oi Cz 


AM_414jj4 


O 1 OA 

z.lzU 


A ni 1^ 
U.Ui jj 


signaling 


Mrrj 


VA Jf A^ Qt f\Q 

AJVl_41oiUo 


1 001 

i.yzi 


U.U4yi 


signaling 


ijKrLMI 


iNM_ZU4y /o 


1 Qn'2 


n noon 
U.UzzU 


signaling 


CZUOIJjZ 


VA/T A 1 nACiCi 


1 QQn 


n n 1 1 o 
U.Uli / 


signaling 


CKAnrl 


MAT r\Ai A'incon 

JN M_uu i u J u D J y 




A A1 1 C 
U.Ullj 


signaling 


LlorjW/ 


MA /I nn 1 r\i i nc i 


V.I l\J 


A (\A 1 1 

U.U411 


signaling 




AM_41 juoj 


1 /^^T 
i.OJ / 


nni nn 
U.UI / / 


signaling 


AKHijiLr IZ 


VA^ /1 1 TO on 
AM_41 /oyU 


1 All 

i.4J / 


A AAI A 

U.UUi4 


signaling 


LI Br J 


VA /f A'^^A A A 

AM_4zd444 


1 /111 


A n-ino 

u.u4uy 


signaling 


rsKjcLt 


IN M_UU 1 U i U o4 i 


1 'Jon 


U.UjU/ 


signaling 


vrfc 1 A 
rUElA 


AM_4ziyDy 




A ni oo 
U.Ulz/ 


signaling 


CKMrSr 


AM_4z46Ui 


1. jjj 


n ni CO 
U.UI jy 


signaling 


QI}/~' A D1 

oKLrArl 


IN M_UU 1 U oU 1 U i 




n CiAO^ 
U.U4ZO 


signaling 




AJVl_41o/4o 


1 1 cm 

i.iyo 


A A/nc 
U.U4/J 


signaling 


ISrrrKZ 


MA* nr\i no a ooc 


1 1 0/l 

i.lz4 


A AOOO 

U.Uzoj 


signaling 


yVlS 14 


iNM_2U4/i5:> 


1 11/1 
1. ii4 


A ni on 
U.UloU 


signaling 


rLrU4 




1 nnn 

i.uyy 


n ni o 1 
U.UlZi 


signaling 


CiUC^Z 


MA>f 'lA/IC/in 

JNM_zU4j4U 


1 noo 


A AOTT 

U.Uj / / 


signaling 


J On I 


MA /I r\r\ 1 nn i /i t 


1 ni c 


U.Uj4/ 


signaling 


KJyUj 


VA4^ All A <Q 

AJVl_4zzi JO 


n QOT 

U.oo / 


n nnn /I 

u.uuy4 


signaling 


AKHUArZo 


VA JI A "i c\i Ar\ 

AM_4iyi4U 


A /IOC 


A AOO 1 

U.Uzjl 


signaling 


Wis ID 


iNM_ZU4D / J 


'X ono 


A f\A OO 

U.U4zZ 


signaling 


LUL,4Zoy0l 


XTA4^ nni 1 AOf^n\ 
JN-IVl_UUli4zu / i 


-J.4Dy 


n nnco 
U.UUjo 


signaling 


VAVZ 


MA 4" OA/1 1 /lO 

JNJVl_zU4i4z 


-3.3y / 


A no 'JO 
U.UZjZ 


signaling 


I rSL-iuzu 


AIV1_UU i Z O J U 1 4 


-J.i j4 


n mcQ 
u.ujoy 


signaling 


III I 


AJVl_4zUo zz 


"2 noA 


n nAA^ 
U.UUoj 


signaling 


rZKiZ 


AM_4z JOD / 


-j.u/y 


A nOQ/l 

u.uzy4 


signaling 


/"'DA DO 

KjKArZ 


AiVl_UU 1 Z j4Uo i 


-z.yo / 


A n/ino 
U.U4UZ 


signaling 


7\TDD/i A 

II\rr4A 


AiVl_4iOooO 


-z. /4y 


n ni 'in 
U.UI jU 


signaling 


1 1 bis 1 


YA4^ /I 1 A7 1 '^ 
AiVl_410 / i J 


0 '7/1/1 
-Z. /44 


n ni 
U.UI jy 


signaling 


iLUlyKrSZ 


IN JV1_ZU4 i ZU 


0 ^i/l/l 
-Z.044 


A n'51 '7 

U.UjI / 


signaling 


rijn 


YA4^ /Ionian 
AM_4zUoDy 


0 /^l 1 


nni 0/1 

u.uiy4 


signaling 




VAJ /lOOOIT 

AiVi_4zyzi / 


0 Al 
-z.olu 


n nn/is 
U.UU40 


signaling 


L,UL-4ZU4Uj 


YAyi A 1 c<no 
Aivi_4 1 o J uy 


0 

-Z.jJi 


A f\r\AQ 

u.uu4y 


signaling 




i> iVl_UU 1 UZ4 J o 3 


0 <;ni 

-Z. Jul 


n non^ 

U.UZUJ 


signaling 


/^DD7Ci 


XTAyf nni non i n< 


-Z.J /O 


n nAQ'7 
u.uuis / 


signaling 




VA* /I 1 "^AAO 

AM_413yyo 


O OCiC 

-Z.JJO 


A AO OA 

U.UZjU 


signaling 


FGF14 


NM 204777 


-2.322 


0.0175 


^iicnalitiO' 


SFRP4 


XM_418831 


-2.204 


0.0162 


signaling 


LOC430333 


XM_001235474 


-1.880 


0.0153 


signaling 


GARNLl 


XM_421244 


-1.850 


0.0477 


signaling 


RAPHI 


XM_421961 


-1.842 


0.0035 


signaling 


C14orfl38 


XM_421460 


-1.756 


0.0012 


signaling 


LOC421876 


XM_419893 


-1.750 


0.0333 


signaling 


NLEl 


XM_415857 


-1.735 


0.0463 


signaling 


CHRM2 


NM_001030765 


-1.721 


0.0310 


signaling 
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Table 3. iii) DE-transcripts related with signaling, 


ion transportation, extracellular matrix protein, and carbohydrate metabolism or post- 


translfltioti P'Ivcos 


vlation moflificatioTi (T^OTitiniied^ 








Gene symbol 


Transcript ID 


Log2 units (strong V S weak) 


p -value 




Category 


LUL/Ooyjo 


AJVl_UU 1 Zjz 1 Zo 


1 ^^99 

- i .oZZ 


U.UZ /o 


signaling 


J? A 1f\Ti 


AiVL_UU 1 ZjjO I J 


-l.j40 


n n/1'7'3 
U.U4/3 


signaling 




AiVl_4 i V J 04 


1 

-1.j4j 


A AA9A 


signaling 


L,UL./Oyjl / 


AiVi_UU 1 ZD 1 y44 


1 ^"51 
-l.jjl 


A A0Q7 
U.UZiS / 


signaling 


rLiAUL-Z 


AiVl_'41oOi J 


-i.jyjj 


A AI 


signaling 




VA/T A 1 /1A99 




A A9'7^ 


signaling 


r^Q PR 


AIVl_UUiZo40Uo 


-1.40 b 


A A/1 QO 


signaling 


LUC4Zyi Oj 


AIV1_4Z0 / i o 


1 A'2.A 
-1 .4:)4 


A A99A 


Signaling 


J3T Y\T/^1 


YA/I Alfi^A'^ 
AiVl_4 101 4 J 


-l.:jyi5 


A AI A'3 


signaling 


KCJM£SU4_I ygy 


YAT A^C^ClQC^ 


-l.o:Sy 


A AA1 Q 

u.uu ly 


signaling 


L>T V\TA 7 


YA4^ /1 1 yl -^nfi 

AJVl_414j /U 


1 'JA^ 

-I.jUj 


A CiA/^l 

U.U40 / 


signaling 


1 jrAiyJ 


YAyl /10n/^</1 
AiVl_4ZU0j4 


-V.zl 1 


A AI AO 
U.U14Z 


signaling 




iN iVl_UU 1 1 Z / i / 1 


-1 .zzo 


A A/1 1 /I 
U.U414 


signaling 


\ \TV A in 


AiVl_UU 1 Z J :>00 1 


-1.1/0 


A A/1/1 1 
U.U441 


signaling 


K\^JMiSUh-_1 oCll 


IN jvi_uu 1 u 1 zyuv 


-i.ijj 


A A/1'5*^ 


signaling 


orKrLUZ 


YA4^ /llQ^/11 
AIV1_4 1 yo4 1 


-1 . iZ4 


A AI /I A 


signaling 


MUdKLIA 


YA4^ /19n^A1 

AM_4ZU0U1 


1 1 1 
-l.llo 


A AA9/i 

U.UUZo 


signaling 


ALiijZL'L, 


YT? nO/^CT^ 
AK_UZOo /O 


1 H'^ft 
-l.U:?o 


A AI 1 0 

u.uiiy 


signaling 


A/fDD 7 

Mrrl 


JNJVi_uuiuu/yi / 


-1.U30 


A A/lO/1 

U.U4Z4 


signaling 


rLrr Iz 




-l.UlO 


A AAA/^ 


signaling 


\M 7 3 


YA/T /I 1 ^AOA 

AM_41 jOoo 


-l.UU/ 


A A 1 A^ 


signaling 


A DT 7/1 


YAyf /11/K<9 

AM_414jjZ 


-yj.y Id 


A A/IQO 

u.U4yz 


signaling 


AUKAZD 


AJVI_4ZjZUj 


(\ Q0< 


A A/I 1 ^ 

U.U41 J 


signaling 


CCKAK 


MAfl" f\n 1 no 1 1 


n o*^o 
-U.oOo 


A A9<9 
V.VJJZ 


signaling 


Ki^J MiSUH-_jfll D 


IN M_UU I vjyjyyjz 


n Q1 < 
-U.ol J 


A A'5'7'3 


signaling 


'T'Or^ 7 T\'^A 

IdL,IL>zh- 


YA4^ r\n 1 9^9 9n/i 

A M_uu 1 z J z zy o 


A A9n 

4.ozy 


A AA9 9 


11 




YA4^ y199r\91 

AJVl_4ZZUZ I 




A AA/19 


TT 
11 


Ki^iyl Z 


AiVl_4Z00 14 


9 n'7'7 


A AA1 < 


11 






9 n9n 


A A9Q9 

v.Vjyz 


TT 
11 


\TTDA T A 


YA4^ A \ A 

AJVI_414jOO 


9 90<^ 

z. /yo 


A AI A9 


TT 
11 


Ai rO VVAh- 


MA4^ nm ncni no 
IN JV1_UU lUoU 1 uz 


z. /yu 


A A'^QT 


TT 
11 


LrKIlyZD 


YAyr /1 1 ^oc\A 
AiVl_4 1 0ZU4 


0 7Q1 

Z. /ol 


A A0Q1 

u.uzyi 


TT 
11 


KLIyJl 


YAT /19C7nC 

AM_4Z J / y J 


Z. jOJ 


A AA1 1 


TT 
11 


rUK 


YA4^ /I 1 ^7/^0 
AiVl_41 J /Oo 


9 9Q/1 

z.zy4 


A A A/1 9 


TT 
11 


Ki^PiKZ 


AiVl_UU 1 z j4zoy 


0 O'yn 
z.z/u 


A A/1 /IT 
U.U44 / 


TT 
11 


KIKKLLj 


YD n9/^0'7/1 

AK_UZOo /4 


0 1/19 


A AI Q9 


TT 
11 


kjAdKkjZ 


iN iVl_ZUD J 4 J 


1 QQ< 


A AI 

U.Ul /o 


TT 
11 


cr r^A A 1 


IN JV1_ZU J J ZZ 


1 TQ/I 

1. /y4 


A AA/1 1 


TT 
11 


iyUUrA/ 


YAyf /1 1 0 1 0^ 
AiVl_4 1 o 1 o J 


u.yyj 


A C\A(^A 
U.U404 


TT 
11 


JrtlD 


YA4^ A^ A^ 00 
AiVl_4 141 y z 


n Qo^ 


A AI O'J 


TT 
11 




AIV1_UU 1 Zj 1 zo J 


r» 8^0 


A A/100 


11 


Jzr CAD J 


YA/T A^ '^S^^^ 
AiVl_4 i J o J J 


-J. jyj 


A CVXAO 


11 


SPATA22 


XM_00 1235 167 


-3.590 


0.0176 


IT 


RCJMB04_lfl 


NM_001031133 


-3.197 


0.0328 


IT 


LOC428404 


XM_425965 


-3.177 


0.0442 


IT 


ATP13A3 


XM_422709 


-3.059 


0.0065 


IT 


GABRB2 


XM_00 1232377 


-2.963 


0.0472 


IT 


SLC8A1 


NM_001079473 


-2.694 


0.0342 


IT 


CACNA2D3 


XM_414338 


-2.622 


0.0100 


IT 


SERINC5 


XM_424762 


-2.615 


0.0129 


IT 
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Table 3. iv) DE-transcripts related with signaling. 


ion transportation, extracellular matrix protein, and carbohydrate metabolism or post- 


translation glycosylation modification (Continued) 








Gene symbol 


Transcript ID 


Log2 units (strong V S weak) 


p -value 


Category 


LUC4ZIO00 


Vr> AT7 ^ AQ 

AK_(Jz / i4o 


-z.ol:S 


A (\1f\l 


TT 
11 




AM_4Z3U/3 


-Z.JiSO 


A A01 0 


TT 
11 


LUL/ /Zjyi 


aJV1_UU i Z o j j J j 


-Z.J / J 


A Al 


11 


KL-jy/KI / 


AJVl_4iy4/ / 


-z.jUd 


A no/n 
U.UU4 / 


TT 
11 


KClIJIo 


VAff /I "ICO IT 

AM_4z!)Zi / 


-z.Juj 


A AA'7|C 

U.UU/O 


TT 
11 


LNNMl 


VA T /111 "Tm 

AM_4zi / Uj 


o "snn 
-z. juy 


A A/l A 

U.U44o 


TT 
11 


KCJyJj 


AJVl_4i /oo4 


-z. 135 


A 

U.U3Z4 


TT 
11 


KCJM£Su4_IleI U 


rsl M_UU 1 UoUoo u 


o 1 no 
-z.iUz 


A no OA 


TT 
11 


ijKlrsoA 


VA Jf C\f\ 1 T2 O 1 O 1 

AM_UU 1 z J z 1 o i 


-i.u<5U 


A AQCC 


TT 
11 


AirO VIA 


MA 4^ '^r\ACin A 


1 ono 


A nAA< 


TT 
11 


D/^ 7\4T>n/1 7 7 "7 

KCJM£Su4_l OUl z 


XTAfl^ (\f\ 1 r\Q 1 Qr\< 
rsl M_UU 1 U J 1 J 


1 n^n 


A noAi 
U.UZol 


TT 
11 


L-jy/LrAJ 


MA/f 0O<001 

JNM_zU3zzl 


-u.yoz 


A C\AlQk 


TT 
11 


DO DV/i 

rZKX4 


rNM_zU4z9i 


-U.KoZ 


A AO^O 


TT 
11 


j| 4 Z7/^ TP T f\ 

MEUr lU 


VA4 A1 Al 1 n 

AJVl_4z4/iy 




A A^OO 


TjA^ 


rAlZ 


VA X A AC^O A 

AM_414jo4 


j.oyZ 


A AA<0 


TjA^ 


^uKZ 


XTA4^ on /I CIO 


z.yuj 


A AAQT 


TjA^ 

liM 


\TT>V\73 

JyKAJyj 


VA4^ A^^1Cl^ 
AM_4z 1 Zy / 


Z.04J 


A A1 

u.ui jy 


TjA^ 


T A TiJf A A 

LAMA4 


vAif Aianof^ 


z.joy 


A A'501 
U.UJZZ 


rLiVl 


\T'T'AT/~' 7 

IS USUI 


VA 4 An 1 T2 1 A A£. 

AJVl_UU 1 Z J 1 44d 


o nn/1 
Z.UU4 


A A'20Q 


TjA^ 




XTA A nn 1 n 0 m 1 1 


1 mn 


A A 1 QO 


TjA^ 

liM 


LAMCl 


MA 4^ On-1 1 /i/i 


l.OoU 


A AO'n'? 

U.UZ J / 


rLiVi 


T A HyfDA 

LAMB4 


VA 4^ nn 1 '1 onn 
AM_UUlzjzo / / 


i.04Z 


A Al Al 

U.Uiul 


TjA^ 


DTHTTA 1 

rrrlAl 


VA4 A'^ir\'7A 

AM_4zlU/4 




A Al AA 
U.UlDO 


TjA^ 


CLVJSzU 


VA/f nniT2onm 
AM_UU 1 ZjZUVZ 


1 A'XQ 


A A 1 n 1 
U.UIUI 


TjA^ 

liM 


L-JJMy 


AM_UU iZji j\j i 


i.zyo 


A A'21 1 


TjA^ 


CM AD 


VA/T A 1 <T2/i 

AM_4iDZoO 


l.zuo 


A AO/f A 


rLiVl 


LULjyOUZO 


XTAfl^ one 1 TO 




A A/im 


rLiVi 


rCutlZl 


JNM_UU1UU1 /jy 


1.1/4 


A At QT 
U.Ulo / 


liiVl 


77 n r"\ D 7 

EruKi 


VAyf A 1 OOQA 


1 1 <c 

l.i JO 


A A/1 An 


TjA^ 


L-rJSho 


VA/T nn 1 1 HOC 
AM_UU i z J 1 3 o o 


1 1 O 1 

i.lzi 


A AO 1 n 

u.uziy 


TjA^ 


/^/~\T A 1 

ajLlZAl 


MA/f oncn'^i 


u.yyo 


A AAA'3 


TnAyf 

rLiVi 


FKrZ 


AM_41oJoz 


A QQ'i 

KJ.yoj 


A C\AA1 


TjA^ 

liM 


CU/Z 


MA 4 'inCACO 


n Q<< 


A AOA'Z 


TjA^ 

liM 


JSIiSJz 


VA4 /I 1 i^lQI 

AJVl_41oJoz 


/I n^n 


A A/inn 


TjA^ 

liM 


UlUr 


VA^ /1'IAAI C 

AM_4zUUi!) 


-j.you 


A AA1 A 


TjA^ 

liM 


CULijAI 


VA /T nn 1 o ^ 1 o 
AM_UU 1 z Jzz 1 o 


Q o^n 
-J.ZoU 


A AAAA 


TjA^ 

liM 


C1/I7D 7 


VA/T /lO /I ni T 

AM_4z4y i / 


-z.yo / 


A AO^A 


liM 


UlUri 


VA^ /lomnn 
AM_4zU /yu 


-Z.ojV 


A AOA'3 
U.UZUJ 


liM 


LrriSiVltS 


VAJ A1'^QQ^ 

Aivi_4Z jyy 1 


-Z. / / 1 


A (Wll 


liiVl 


PVP J 
rKrl 


YA4 /iiQi/in 
AiVl_4 1 y Z4U 


-z.j'+y 




liiVl 


UxiviAZ 


VAyr A-iQlAf^ 
AiVl_4 1 y / 40 


-Z.jjy) 


A AAAA 


rLM 


COL8A2 


XM 425780 


-z.zyj 


A ni 7n 


rLM 


FNBP4 


XM_424260 


-2.279 


0.0486 


EM 


EGFL6 


XM_416835 


-2.268 


0.0392 


EM 


CDH18 


XM_426046 


-1.955 


0.0090 


EM 


CLDN8 


XM_425544 


-1.598 


0.0288 


EM 


RCJMB04_34k20 


NM_001031112 


-1.214 


0.0085 


EM 


SRPX 


XM_416781 


-1.123 


0.0244 


EM 


DLGl 


XM_422701 


-1.084 


0.0171 


EM 


FBLNl 


NM_204165 


-0.981 


0.0206 


EM 
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Table 3. v) DE-transcripts related with signaling, ion transportation, extracellular matrix protein, and carbohydrate metabolism or post- 
translation glycosylation modification (Continued) 



Gene symbol Transcript ID Log2 units (strong VS weak) p-value Category 



F13A1 


JNM_zU4doj 


-0.949 


0.0234 


EM 


FREMl 


XM_424932 


-0.910 


0.0057 


EM 


MhCjriU 


XJVl_4Z4/iy 


3.910 


U.0398 


EM 


MGAT4C 


XM_425447 


3.007 


0.0497 


/^A if" /~1A if 

OM or CM 


CHST3 


NM_205121 


2.839 


0.0372 


/~1A if /~1A if 

GM or CM 


tUtMJ 


AJVI_4ZZ293 


2.539 


0.0179 


UM or CM 


LARGE 


JN M_UU i 0U44U4 


2.071 


0.0429 


/~'A>f /"'AT 

LrM or CM 


GtFlZ 


aM_4245 /i 


1.919 


A A1 '71 

0.0173 


/^A >ff /^A >f 

(jM or CM 


/—I A T K TT'T 1 


AM_UU 1231 964 


1 one 

1.895 


0.0452 


/^A >ff /^A >f 

CjM or CM 


WDR77 


NM_00 10309 16 


1.805 


0.0445 


/^A if /^A if 

GM or CM 


NDST3 


"\7"A if /I y"> C 

XM_426325 


1.403 


0.0121 


/~"A if /~1A if 

GM or CM 


BJGALl 1 


"VAyl Ai^^COA 

aM_42ojo4 


1.254 


A A/( 01 

0.0483 


y^A if y^A >f 

CjM or CM 


OGDHL 


VAT /n 1 crxT 
aM_4201)3 


1.144 


A AA-no 

0.0022 


uM or CM 


MANIAZ 


VA J Ai£.Ar\n 
AJVl_41o4yU 


-4.700 


A AAC? 

0.0087 


/^AT /^A T 

CM or CM 


POFUT2 


XM_421892 


-2.804 


0.0156 


/~1A if /~1A if 

GM or CM 


LUL/ /Z1d4 


aJVI_UU i Ijj 32y 


-2.798 


0.0006 


/^A/T /"'AT 

uM or CM 


KLB 


aM_423224 


-2.335 


A AO OA 

0.0380 


CM or CM 


TRIM/. 2 


NM_001099354 


-1.986 


0.0041 


/~1A if /~1A if 

GM or CM 


LOL/ /iJOi 


AJVl_UU 1 23464 / 


-1.433 


0.0369 


CM or CM 




i> LVL_\J\J 1 V J 7 J i U 


-1.324 






NDST4 


XM_420638 


-1.239 


0.0261 


GM or CM 


PFKM 


NM_204223 


-1.170 


0.0460 


GM or CM 


NUP153 


XM_4 18937 


-1.109 


0.0378 


GM or CM 


GPDIL 


XM_4 18763 


-0.947 


0.0383 


GM or CM 


PHKA2 


XM_416811 


-0.926 


0.0326 


GM or CM 


B3GNTL1 


XM_415599 


-0.883 


0.0109 


GM or CM 


PMMl 


XM_416228 


-0.787 


0.0490 


GM or CM 


MMPll 


XM_00 1232776 


2.209 


0.0391 


GM or CM 


ST3GAU 


XM_417860 


-1.043 


0.0094 


GM or CM 



rr represents ion/proton transporter, EM represents extracellular matrix, GM represents post-translation glycosylation modification, and CM represents 
carbohydrate metabolism. 



As shown in Table 4 and Table 5, many DE-transcripts 
are involved in signal transduction, such as GO terms 
purinergic nucleotide receptor activity (GO:0001614), 
nucleotide receptor activity (GO:0016502), purinergic 
receptor activity (GO:0035586), transmembrane signaling 
receptor activity (GO:0004888), and negative regulation of 
BMP signaling pathway (GO:0030514). Among them, 
GO:0004888 dominantly contains 33 transcripts encoding 
signal receptors, and these receptors could be further 
classified into several subgroups: OXTR, LOC431251 and 
SSTR3 belong to reproductive hormone receptors; CHRM2, 
ADRA2B, P2RX4, P2RY2, EDNRB2, GABRB2, GABRG2, 
LOC428961 and NPFFR2 function as receptors mediating 
neurotransmitters or neuropeptide; GRIN2B and GRIN3A 
could modulate the efficiency of synaptic transmission; 
NTRKl and NTRK2 belong to the receptor tyrosine kinase 
(RTK) family, and are involved with neurotrophin 
(GO:0005030 - neurotrophin receptor activity; and 
GO:0043121 - neurotrophin binding) (Table 5). 



Besides various enriched molecular function shown 
above, many biophysical processes are also found to be 
enriched among the DE-transcripts, including a series of 
processes and subgroups (Tables 4 and 5). GO:0003951 
(NAD"^ kinase activity) modulate the metabolism or redox 
in cell (Table 5). Enrichment of GO:0009409 (response to 
cold) may reflect the fact the rearing condition of 
experimental hens was in the winter at room temperature 
about 2 to 10°C. GO:0046209 (nitric oxide metabohc 
process) may regulate vascular or smooth muscle relaxation 
or other functions. GO:0002028 is involved in ion 
transportation. While the subgroup processes of muscular 
development and activity include skeletal muscle fiber 
development (enrichments of GO:0048741, GO:0048747 
and GO:0055002) and striated muscle contraction 
regulation (enrichments of GO:0055117 and GO:0006942). 
It is of note that there is almost no striated muscle in avian 
uterus except smooth muscle. However, the chicken 
genome project was completed in 2004, and the functional 
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Table 4. Enriched gene ontology (GO) terms revealed from identified DE-transcripts according to biological_process ontology 



Group 


GOID 


Term 


P 


Gene symbol or representative pubUc ID 


Reproductive hormone 


GO:0060126 


Somatotropin secreting cell 


0.013 


0TX2, WNT4 


S5Tithesis and regulation 




differentiation 








GO:0021984 


Adenohjrpophysis 


0.044 


0TX2, WNT4 






development 








GO:0032355 


Response to estradiol stimulus 


0.044 


S0CS2, AREGB 


Signal transduction 


GO:0030514 


Negative regulation of BMP 


0.030 


TOBl, GREMl 






signaling pathway 






Biophysical processes 


GO:0048741 


Skeletal muscle fiber 


0.012 


SLC23A2, CHAT 






development 








GO:0015074 


DNA integration 


0.012 


LOC770294, LOC770705, ENS-3 




GO:0055117 


Regulation of cardiac muscle 


0.013 


P2RX4, NKX2-5 






contraction 








GO:0009409 


Response to cold 


0.018 


IL4, SLC27A1 




GO:0048747 


Muscle fiber development 


0.022 


SLC23A2, CHAT 




GO:0046209 


Nitric oxide metaboUc 


0.024 


P2RX4, CPSl 






process 








GO:0007586 


Digestion 


0.030 


PGA5, PRSS2, LOC396365 




GO:0015849 


Organic acid transport 


0.033 


SLC23A2, 0CA2, LOC770309, 










SLC7A14, SLC27A1 




GO:0046942 


Carboxylic acid transport 


0.033 


SLC23A2, 0CA2, LOC770309, 










SLC7A14, SLC27A1 




GO:0055002 


Striated muscle cell 


0.034 


SLC23A2, CHAT, TTN, NKX2-5 






development 








GO:0006942 


Regulation of striated muscle 


0.037 


P2RX4, NKX2-5 






contraction 








GO:0002028 


Regulation of sodium ion 


0.044 


NKX2-5, NEDD4L 






transport 






Reproductive biophysical 


GO:0060748 


Tertiary branching involved 


0.009 


WNT4, AR 


processes 




in mammary gland duct 










morphogenesis 








GO:0060745 


Mammary gland branching 


0.013 


WNT4, AR 






involved in pregnancy 








GO:0060562 


EpitheUal tube 


0.019 


DEAFl, WNT3, GREMl, WNT4, 






morphogenesis 




NKX2-5, HOXAIL AR, AREGB 




GO:0060444 


Branching involved in 


0.020 


WNT4, AR, AREGB 






mammary gland duct 










morphogenesis 








GO:0009994 


Oocyte differentiation 


0.024 


WNT4, GDF9 




GO:0048599 


Oocyte development 


0.024 


WNT4, GDF9 




GO:0060603 


Mammary gland duct 


0.033 


WNT4, AR, AREGB 






morphogenesis 








GO:0060135 


Maternal process involved in 


0.037 


WNT4, AR 






female pregnancy 







GOID represents the identifiers, and Term represents term definitions for Gene Ontology term entities, p: p-value of significance (Welch t-test). 



gene database of G. gallus remains incomplete, some 
ontology annotations of DE-genes may refer to mammalian 
homologs, which may account for our results. The genes 
related to muscular cell contraction are likely to modulate 
the mobility of uterus to facilitate egg rotation and 
calcification (Johnson, 1986; Jonchere et al., 2010). 
Similarly, there is no digestion process in the uterus, three 
genes in GO:0007586 (digestion process), PGA5 (an 



aspartic acid protease, which is involved in ovulation 
(Peluffo et al., 2011)), PRSS2 and LOC396365 
(preprogastrin), are likely to promote the maturation of 
secretary extracellular proteins or regulate the secretion of 
uterus glands and mobility of uterus. 

The final group of reproductive biophysical processes 
also includes several subgroups of processes (Table 4). 
Epithelial tube morphogenesis (GO:0060562) may regulate 
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Table 5. Enriched gene 


ontology (GO) terms revealed from identified DE-transcripts according to molecular function ontology 


Group 


GOID 


Term 


P 


Gene symbol or representative public ID 


Signal transduction 


GO:0005030 


Neurotrophin receptor 
activity 


0.013 


NTRKl, NTRK2 




GO:0001614 


Purinergic nucleotide 

receptor activity 


0.017 


P2RX4, P2RY2, ENSGALG00000012080 




GO:0016502 


Nucleotide receptor 
activity 


0.017 


P2RX4, P2RY2, ENSGALG00000012080 




GO:0043121 


Neurotrophin binding 


0.024 


NTRKl, NTRK2 




GO:0035586 


Purinergic receptor 
activity 


0.026 


P2RX4, P2RY2, ENSGALG00000012080 




GO:0004888 


Transmembrane signaling 
receptor activity 


0.049 


OXTR, LOC431251, SSTR3, CHRM2, 
ADRA2B, P2RX4, P2RY2, EDNRB2, 
GABRB2, GABRG2, LOC428961, 
NPFFR2, GR1N2B, GR1N3A, NTRKl, 
NTRK2, EPHB6, DDR2, TMPRSS6, 
PCSK5, CCKAR, IFNAR2, CSFIR, 
TLR5, ORlOAl, LOC768958, 
LOC769317, LOC777484, GPR39, 
GPR97, ENSGALG00000017405, 
ENSGALGOOOOOOl 7093, 
ENSGALG00000012080. 


Biophysical processes 


GO:0003951 


NAD+ kinase activity 


0.013 


C5orf33, NADK 




GO:0005319 


Lipid transporter activity 


0.049 


ATP lie, ATP8A2, ATP8B3, APOB, 
LOC769564, SLC27A1 



GOID represents the identifiers, and Term represents term definitions for gene ontology term entities, p: p-value of significance (Welch t-test). 



the development of uterus glands (tubular epithelial glands). 
Oocyte development subgroup contains oocyte 
differentiation (GO:0009994) and oocyte development 
(GO:0048599). Female pregnancy subgroup contains 
enrichments of GO:0060135 , 00:0060745 , 00:0060748, 
00:0060444 and 00:0060603. 

Overall, laying is an avian reproductive behavior, and 
eggshell calcification is regulated by relative reproductive 
hormones and neurotransmitters, which may finally affect 
eggshell quaUty through a complex suite of biophysical 
reactions. 

Confirmation of DE-transcripts by qRT-PCR 

21 DE-transcripts (9 up-regulated and 12 down- 
regulated) were chosen for validation using qRT-RCR 
experiments, and the four microarray samples were tested in 
pairs for #19 vs #35 and #40 vs #80, respectively. As shown 
in Figure 3, 16 out of the 21 tested transcripts (76%) were 
confirmed by qRT-PCR experiments, though the absolute 
fold-change values are slightly different. The remaining 
transcripts, CRYBBl, EX0C6B, LOC416916, MAN1A2 and 
CHST3, showed inconsistent differential expression 
between qRT-PCR and microarray experiments. 

GAPDH, CHST3, GALNTLl, NDST4, LARGE, SPl, 
RH0BTB2, and WDR72 were selected to examine the PCR 
efficiency. The results showed the PCR efficiency of these 
genes ranged from 86.8% to 94.2%, and the PCR efficiency 
of inner reference (GAPDH) and other genes seemed nearly 



similar. 

DISCUSSION 

Laying is regarded as avian reproductive behavior, 
which is regulated by reproductive hormones and 
neurotransmitters. The chicken oviduct has been 
extensively used as a model to study hormonal induction of 
protein synthesis (Khuong and Jeong. 2011). Under the 
control of steroid hormones or neurotransmitters, the 
tubular gland epithelial cells synthesize and secrete a great 
variety of proteins to form egg white and eggshell when egg 
passes through the oviduct (Mann et al., 2006). In this paper, 
1,195 DE-transcripts have been identified to be related with 
eggshell strength. OOEAST analysis further identify some 
significantly enriched GO terms, and the enriched GO terms 
suggest that some DE-transcripts mediate reproductive 
hormones or neurotransmitters to affect eggshell quaUty 
(Tables 4 and 5). 

Both terms 00:0060126 and 00:0021984 are involved 
in reproductive hormone regulation, and share two genes, 
0TX2 and WNT4. 

Otx2 is a paired-like homeodomain transcription factor, 
which can mediate GnRH (gonadotropin releasing 
hormone) signaling (Kelley et al., 2000). Functional studies 
revealed that Otx2 is required as early as gastrulation for 
neural induction, and even for brain development (Rhinn et 
al., 1998). However, Otx2 is also of importance for 
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Figure 3. Real-time RT-PCR validation of microarray data. Expression levels of the first 12 transcripts (AMDHDl , ATP6V1A, 
CRYBBl, FOB, GAS2L3, GIT2, MAN1A2, NDST4, PLCXDl, RCJMB04-34k20, SLC8A1 and TBXASl) were down-regulated in 
microarray experiment, and the last 9 transcripts (ACYP2, CA5B, CHST3, C0L12A1, CRABPl, EX0C6B, LOC416916, NPY and 
WDR72) were up-regulated in microarray experiment. Among the above transcripts, five transcripts (CRYBBl, MAN1A2, CHST3, 
EX0C6B and LOC416916) failed to be verified by real-time RT-PCR. 



neurogenesis and cellular proliferation in multiple other 
tissues (Layman et al., 2011). 

As a member of the WNT family, Wnt4 is a secreted 
glycoprotein signaling molecule and involved in paracrine 
signaling (Diaz et al., 2011). Wnt4 is critical for female sex 
determination and differentiation (Chen et al., 2011). In the 
female, Wnt4 is positively involved in ovarian 
development; while in the male mutated WNT4 will result 
in aberrant testis development (Diaz et al., 2011; 
Barrionuevo et al., 2012). On the other hand, Wnt4 is also 
potent to regulate the development of the female 
reproductive tract (Franco et al., 2011). Furthermore, WNT4 
is expressed postnatally in ovarian follicles and corpora 
lutea, and its expression increases in response to 
gonadotropin (Hsieh et al., 2002). Wnt4 mediates follicle 
development and fertility by regulating the expression of 
genes involved in steroidogenesis, prostaglandin 
biosynthesis, tissue remodeling, and angiogenesis (Hsieh et 
al., 2002; Boyer et al., 2010). Moreover, excluding its 
reproductive contributions, WNT4 is also tightly associated 
with bone strength (Zmuda et al., 2011). 

Our results also show that some DE-transcripts are 
involved in signal transduction (Tables 4 and 5), among 
which, NTRKl, NTRK2, P2RX4, and P2RY2 are overlapped 
in multiple enriched GO terms (Table 5). 

Ntrkl, also named TrkA, and Ntrk2 TrkB, are two 
members of the neurotrophic tyrosine kinase receptor 
(NTKR) family. These kinases are membrane -bound 



receptors mediating various functions of neurotrophins, 
such as cell survival, migration, outgrowth of axons and 
dendrites, synaptogenesis, remodeling of synapses, and 
synaptic transmission (Ohiral and Hayashi, 2009). So far, 
several neurotrophins have been well studied, such as nerve 
growth factor (NGF), brain-derived neurotrophic factor 
(BDNF), neurotrophin 3 (NT-3), and NT-4/5. NTKRs are 
high affinity receptors of neurotrophins. TrkA mediates the 
biological response of NGF, while BDNF and NT-4/5 are 
the preferred ligands for TrkB (Ohiral and Hayashi, 2009). 
Additionally, NTKRs also play roles in some biomaterials. 
NT-4 may modulate proliferation and differentiation of the 
dental epithelium and promote production of the enamel 
matrix via the TrkB-MAPK pathway (Yoshizaki et al., 
2008). 

Both P2RX4 and P2RY2 are purinergic receptors. 
Purinergic receptors are subdivided into metabotropic P2Y 
receptors and iono tropic P2X receptors. P2Y receptors are 
coupled to G-protein and trigger inositol 1,4,5 -triphosphate 
(IP3)-induced intracellular Ca^^ release following activation 
of phospholipase C, while P2X receptors are ligand-gated 
ion channels. P2RX4 will be discussed later, while P2RY2 
performs a dominant role in calcium signaling during 
osteoblast differentiation (Nishii et al., 2009). It is known 
that extracellular ATP, UTP, and PPi can strongly block the 
mineralization of matrix nodules, while this potent 
inhibition of bone formation is mediated by P2RY2 (OiTiss 
et al., 2007). Furthermore, P2RY2 is also involved in 
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inhibition of intercellular communication between 
osteoblasts (Hoebertz et al., 2003). 

At present, there are at least three cDNA microarray 
studies globally investigating the gene expression in 
chicken shell gland (Yang et al., 2007; Dunn et al., 2009; 
Jonchere et al., 2010), but the overlap among the DE-genes 
from these studies is not plentiful. Different animals, tissue 
samples or treatment methods may partially account for this 
problem. 

Yang et al. (2007) harvested uterus tissues at 2 h post 
oviposition, and screened out 34 known genes in the shell 
glands between hens with low and high egg productive rates. 
This study and our data share a single gene CALDl 
(caldesmon 1) (Figure 4). CALDl is a ubiquitous actin and 
calmodulin binding protein, and functions as a substrate for 
mitogen-activated protein kinase (Childs et al., 1992) or as 
serine and threonine kinases (Sutherland et al., 1994). 

Dunn et al. (2009) identified 266 DE-genes in shell 
glands from 25-week old mature hens comparing to 12- 
week old juveniles from high and low bone quality lines, 
respectively. The tissues sampled when eggs passed 
through the oviducts but not in shell glands. Three DE- 
genes are also found in our data: NADK (NAD kinase), 
LOC422993 (Similar to interferon-induced membrane 
protein Leu- 13/9-27), and LAMPS (lysosomal-associated 
membrane protein 3) (Figure 4), suggesting potential crucial 
function of these genes in not only early stage of eggshell 
calcification but also other stages of eggshell formation. 

Jonchere et al. (2010) used the 40- week old hens at 18 h 



post oviposition (corresponding to the rapid phase of 
calcification), and identified 469 DE-known genes in uterus 
versus both white isthmus, and magnum. There are 7 genes 
consistently identified in their study and our data, such as 
P2RX4 (purinergic receptor P2X, ligand-gated ion channel, 
4), FSTLl (foUistatin-hke 1), TUBGCP4 (Tubulin, gamma 
complex associated protein 4), WDR77 (WD repeat domain 
77), RCJMB04_6gl6 (microtubule-associated protein 1 
light chain 3 beta), PWPl (PWPl homolog in S. cerevisiae) 
and SGKl (serum/glucocorticoid regulated kinase 1) 
(Figure 4). 

On the other hand, three additional DE-genes in our data 
were previously found in the acid soluble part of chicken 
eggshell organic matrix (Mann et al., 2006). These three 
genes, FSTLl (follistatin-like 1), CAMK2D (calcium/ 
calmodulin-dependent protein kinase (CaM kinase) II delta) 
and KRT75 (keratin 75) (Figure 4), could reflect potential 
interaction of eggshell calcification and organic matrix 
formation. 

Among these overlapping DE-genes, both P2RX4 
(Jonchere et al., 2010) and NADK (Yang et al., 2007) are 
also present in our enriched GO terms (Tables 4 and 5), and 
FSTLl (Mann et al., 2006; Jonchere et al., 2010) occurs in 
more than three relative studies. 

P2RX4 is one member of the P2X receptors (P2RX). 
P2RX are ionotropic ATP-gated ion channels conducting 
Ca^"^ inflow (Fodor et al., 2009), with high capability of 
Ca^"^ permeabilities corresponding to at least 100-fold those 
of Na"^ (Burnashev, 1998). In chondrogenic mesenchymal 
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Figure 4. MA-plot of all microarray tested chicken-genes. X-axis: the average normalized expression values across all 4 (strong+weak) 
eggshell samples (in Log2 scale); Y-axis: the log-odds ratio between the average expression values of strong vs. weak eggshell samples; 
grey dots: nondifferentially expressed genes (non-DE genes); red dots: differentially expressed genes (DE-genes) identified in this study; 
big blue dots: DE-genes (CALDl, NADK, LOC422993, LAMPS, P2RX4, FSTLl, TUBGCP4, WDR77, RCJMB04_6gl6, PWPl and 
SGKl) reported in previous studies. 
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cells, P2X4 receptors could conduct Ca^^ inflow to elevate 
intracellular Ca^* levels, and finally promoting extracellular 
matrix production (Fodor et al., 2009). Eggshell 
calcification requires considerable ion transportation, 
especially Ca^"^, and various extracellular matrix synthesis 
and secretion, whether and how P2RX4 channels regulate 
these processes requires further studies. 

NAD kinases (NADKs) are a family of enzymes 
transferring a phosphate group from ATP to NAD to 
generate and maintain the cellular NADP pool (Pollak et al., 
2007). It is reported, that during development of placenta, 
the expression level of NADK appears drastically elevated 
(Lerner et al., 2001). 

Fstll is a secreted glycoprotein belonging to the BM- 
40/SPARC/osteonectin family containing both calcium- 
binding domain and Follistatin-like domain (Hambrock et 
al., 2004). As a mesenchymal factor, Fstll is critical for 
oviduct development, and determines the differentiation of 
secretary epithelial cells and ciliated epithelial cells in the 
oviduct (Umezu et al., 2010). This means Fstll may 
modulate chicken endometrium development during 
eggshell formation. However, Fstll is also present in the 
organic part of eggshells (Mann et al., 2006), and Jonchere 
et al. (2010) propose it may be a uterine antiprotease 

IMPLICATIONS 

Above all, using Affymetrix Chicken Array, 1,195 DE- 
transcripts were identified in the shell gland between "high 
shell strength" and "low shell strength" hens, which 
represent 941 unique known genes. According to gene 
ontology annotations, these transcripts are involved in a 
wide range of biological processes; the most prominent DE- 
transcripts relate to signal transduction, metabolism, 
extracellular matrix, or ion transport and homeostasis, and 
so on. Furthermore, Gene Ontology (GO) term enrichment 
of DE-transcripts suggests that avian eggshell calcification 
is Ukely to be regulated by relative reproductive hormones 
and neurotransmitters, which may finally affect eggshell 
quality through a complex suite of biophysical processes. 
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